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DETERMINING THE FRACTAL DIMENSION 
OF THE INTERSTELLAR MEDIUM 
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RESUMEN 

El Medio Interestelar parece tener una estructura fractal subyacente, la cual se puede caracterizar a traves de 
su dimension fractal (Df). Sin embargo, diversos factores pueden afectar la determination de Df, tales como 
distorsiones causadas por la proyeccion, baja resolution de la imagen, la opacidad de la propia nube, y una baja 
razon senal/ruido. En este trabajo utilizamos tanto nubes simuladas como mapas reales de nubes moleculares 
para estudiar estos efectos y poder estimar de una manera liable la dimension fractal. Los resultados indican 
de una forma autoconsistente que la dimension fractal del Medio Interestelar esta en el rango 2.6 < Df < 2.8, 
lo cual es significativamente mas alto que el valor Df ~ 2.3 que usualmente se asume en la literatura. 

ABSTRACT 

The Interstellar Medium seems to have an underlying fractal structure, which can be characterized through 
its fractal dimension (Df). However, several factors may affect the determination of Df, such as distortions 
due to projection, low image resolution, opacity of the cloud, and low signal-to-noise ratio. Here we use both 
simulated clouds and real molecular cloud maps to study these effects in order to estimate Df in a reliable 
way. Our results indicate in a self-consistent way that the fractal dimension of the Interstellar Medium is in the 
range 2.6 < Df < 2.8, which is significantly higher than the value Df ~ 2.3 usually assumed in the literature. 
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1. INTRODUCTION 

The fractal dimension (Df) seems to be a sim- 
ple but useful tool for characterizing the degree of 
complexity of the Interstellar Medium (ISM). Most 
of studies use the perimeter-area method to calcu- 
late the fractal dimension of the boundaries of pro- 
jected molecular clouds (D per ). The observational 
evidence can be summarized by saying that D per has 
a more or less universal value around ~ 1.35 (see 



Sanchez et al.l 120051 and references therein) . From 



this, it is usually assumed that the real 3-dimensional 
fractal dime nsion is Df = D per + 1 ~ 2.35 (e.g., 



BeechJ 11992J) , a value which moreover is consistent 
with turbulent diffusion in an incompressible fluid 
( Meneveau fc Sreenivasan ll99Ch . In this contribu- 
tion we show that this simple relation between Df 
and D per is not valid. We study carefully the effects 
that different factors have on the determination of 
Df in order to calculate in a highly reliable way its 
value for the ISM. 

2. METHOD 

Here we follow an empirical approach. We first 
simulate 3-dimensional fractal clouds with a given 
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w ell-defined fractal dimen sion Df (see details in 
inchez et al. 2005. 2006). Then we project the 
clouds to generate 2-dimensional maps. When do- 
ing this, we consider and analyze in detail the ef- 
fects produced by the projection itself, image reso- 
lution, cloud opacity, and low signal-to-noise ratio 
(S/N). We take all these factors into consideration 
to develop algorithms to estimate fractal dimensions 
in a precise and accurate way. We use two differ- 
ent fractal estimators to corroborate the results and 
trends obtained: the perimeter-area-based dimen- 
sion (D per ) and the mass-size dimension (D m ). 

3. FRACTAL DIMENSION ESTIMATION 
Our main results can be summarized as follows. 

• Projection. The projected perimeter dimen- 
sion D per decreases as the 3-dimensional fractal di- 
mension Df increases. This can be seen in Figure Q] 
for fractal clouds having dimensions in the range 
2.0 < Df < 2.9. The interesting point here is that 
the observed average value D per — 1-35 is more con- 
sistent with Df~ 2.6 than with Dt ~ 2.3. 

• Resolution. We have decreased the resolution 
by hand and we have obtained that D per decreases 
as resolution decreases. This is an expected result 
since as the pixel size is larger (worse resolution) the 
perimeter becomes smoother because the fractal fea- 
tures (the irregularities in the cloud contours) blend 
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Fig. 1. Perimeter dimension D per as a function of the 
fractal dimension Df calculated for 50 random clouds. 
Each data point is the average value and the bars are the 
standard deviations. Fhe horizontal solid line indicates 
the value D per = 1.35. 

with each other. The results shown in Figure Q] cor- 
respond to "high resolution" images, i.e. cloud for 
which the size is N p i X > 400 pixels. In this case 
we get D per — 1-36 ± 0.03 for the case Df = 2.6, 
but if the resolution is decreased to N p i X = 200 then 
D per = 1.28 ± 0.02. 

• Opacity. We have generated cloud maps with 
different values for the total optical depth r. We 
have obtained the important result that opacity does 
not affect the estimation of D per . This is because 
even for optically thick clouds, for which central and 
densest cores are difficult to detect, the fractal prop- 
erties of the cloud are contained in the shape of the 
external contours. The situation is different for the 
mass dimension D m because this estimator use in- 
formation from all the cloud structure (mass vs. ra- 
dius). This method (D m ) fails when r > 1 but in 
any case it can be used without problems for opti- 
cally thin regions. 

• Noise. Very noisy maps can seriously affect the 
estimation of both D per and D m , artificially increas- 
ing the structure irregularities and therefore decreas- 
ing the final value of Df. This effect can be mini- 
mized by smoothing the map before calculating Df. 
We have shown that this previous step should be 
done in low-S/N maps only if the smoothing process 
maximizes S/N. 

4. APPLICATION TO THE ISM 

Taking all the above effects into account we 
have calculated the fractal dimension of Ophiuchus, 
Perseus, and Orion molecular clouds using emission 
maps in different lines. All the results are summa- 
rized in Figure [2j where we can see that the fractal 



Fig. 2. Perimeter dimension D per obtained for each 
molecular cloud map labeled in the horizontal axis. The 
horizontal lines indicate values calculated for simulated 
clouds with Df values from 2.0 to 2.9 in increments of 
0.1. The interval 2.6 < Df < 2.8 is emphasized with 
thick solid horizontal lines. 

dimension is always in the range 2.6 < Df < 2.8. The 
mass dimension met hod yielded very similar results 
(|Sanchez et alj|2007l) . 

These results are in perfect agreement with re- 
cent simula tions of compressively driven turbulence 
in the ISM ( Federrath et al.ll2007 ). and this is an im- 
portant point because, as pointed out by these au- 
thors, the ISM i s far t o be an incompressible fluid. 
Federrath et all (|2007l) obtained Df ~ 2.6 - 2.7 in 
their standard simulations and they showed that 
Df can be as small as ~ 2.4 only for the extreme 
case of high turbulent Mach numbers (of the or- 
der of 10). Another important point is that it has 
been argued that a fractal ISM with Df ~ 2.3 
could account for some observed c loud properties 
(e.g., lElmegreen fc Falgarond Il996f ). but we have 
also shown that Df ~ 2.6 is roughly consistent with 
the average properties of th e ISM, in particular the 
mass and size distributions (jSanchez et al]|2006h . In 
summary, it seems that Df ~ 2.7 ± 0.1 for the ISM. 
The relevance of this relatively high fractal dimen- 
sion value has to be analyzed, mainly concerning the 
physical processes involved in the structure of the 
ISM. 
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